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ABSTRACT 

This work presents a novel Initial Ranging scheme for orthog- 
onal frequency-division multiple-access networks. Users that 
intend to establish a communication link with the base sta- 
tion (BS) are normally misaligned both in time and frequency 
and the goal is to jointly estimate their timing errors and car- 
rier frequency offsets with respect to the BS local references. 
This is accomplished with affordable complexity by resorting 
to the ESPRIT algorithm. Computer simulations are used to 
assess the effectiveness of the proposed solution and to make 
comparisons with existing alternatives. 

1. INTRODUCTION 

A major impairment in orthogonal frequency-division multiple- 
access (OFDMA) networks is the remarkable sensitivity to 
timing errors and carrier frequency offsets (CFOs) between 
the uplink signals and the base station (BS) local references. 
For this reason, the IEEE 802. 16e-2005 standard for OFDMA- 
based wireless metropolitan area networks (WMANs) speci- 
fies a synchronization procedure called Initial Ranging (IR) in 
which subscriber stations that intend to establish a link with 
the BS transmit pilot symbols on dedicated subcarriers using 
specific ranging codes. Once the BS has detected the pres- 
ence of these pilots, it has to estimate some fundamental pa- 
rameters of ranging subscriber stations (RSSs) such as timing 
errors, CFOs and power levels. 

Initial synchronization and power control in OFDMA was 
originally discussed in [l] and Q while similar solutions can 
be found in lO-H. A different IR approach has recently been 
proposed in |l5l|. Here, each RSS transmits pilot streams over 
adjacent OFDMA blocks using orthogonal spreading codes. 
As long as channel variations are negligible over the ranging 
period, signals of different RSSs can be easily separated at 
the BS as they remain orthogonal after propagating through 
the channel. Timing information is eventually acquired in an 
iterative fashion by exploiting the autocorrelation properties 
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of the received samples induced by the use of the cyclic prefix 
(CP). 

All the aforementioned schemes are derived under the as- 
sumption of perfect frequency alignment between the received 
signals and the BS local reference. However, the occurrence 
of residual CFOs results into a loss of orthogonality among 
ranging codes and may compromise the estimation accuracy 
and detection capability of the IR process. Motivated by the 
above problem, in the present work we propose a novel rang- 
ing scheme for OFDMA networks with increased robustness 
against frequency errors and lower computational complexity 
than the method in |5|. To cope with the large number of pa- 
rameters to be recovered, we adopt a three-step procedure. In 
the first step the number of active codes is estimated by re- 
sorting to the minimum description length (MDL) principle 
|6|. Then, the ESPRIT (Estimation of Signal Parameters by 
Rotational Invariance Techniques) f7 ] algorithm is employed 
in the second and third steps to detect which codes are ac- 
tually active and determine their corresponding timing errors 
and CFOs. 



2. SYSTEM DESCRIPTION AND SIGNAL MODEL 

2.1. System description 

We consider an OFDMA system employing N subcarriers 
with index set {0, 1, . . . , TV — 1}. As in [5], we assume that a 
ranging time-slot is composed by M consecutive OFDMA 
blocks where the N available subcarriers are grouped into 
ranging subchannels and data subchannels. The former are 
used by the active RSSs to complete their ranging processes, 
while the latter are assigned to data subscriber stations (DSSs) 
for data transmission. We denote by R the number of ranging 
subchannels and assume that each of them is divided into Q 
subbands. A given subband is composed of a set of V adja- 
cent subcarriers which is called a tile. The subcarrier indices 
of the gth tile (q = 0, 1, . . . , Q — 1) in the rth subchannel (r = 
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0,1,..., R~ 1) are collected into a set J'q 



(r) 



where the tile index iq ' can be chosen adaptively according 
to the actual channel conditions. The only constraint in the 



(r) 

selection of iq ' is that different tiles must be disjoint, i.e., 

Jq^'^ n jj;'^ = for qi 7^ (72 or ri 7^ ra. The rth ranging 
subchannel is thus composed of QV^ subcarriers with indices 
in the set J^'') = U^Jp^ J'J''^ while a total of Nr = QVR 
ranging subcarriers is available in each OFDMA block. 

We assume that each subchannel can be accessed by a 
maximum number of K^na^ = min{y, M} — 1 RSSs, which 
are separated by means of orthogonal codes in both the time 
and frequency domains. The codes are selected in a pseudo- 
random fashion from a predefined set {Co, Ci, . . . , Ck„^^-i} 
with 

[C,].,™ = e^-2-'=(v-T+^) (1) 

where v — 0,1,...,^ — 1 counts the subcarriers within a 
tile and is used to perform spreading in the frequency do- 
main, while TO = 0, 1, . . . , Af — 1 is the block index by which 
spreading is done in the time domain across the ranging time- 
slot. As in 0, we assume that different RSSs select different 
codes. Also, we assume that a selected code is employed by 
the corresponding RSS over all tiles in the considered sub- 
channel. Without loss of generality, we concentrate on the rth 
subchannel and denote by K < Kn-mx the number of simulta- 
neously active RSSs. To simplify the notation, the subchannel 
index *^'''' is dropped henceforth. 

The signal transmitted by the fcth RSS propagates through 
a multipath channel characterized by a channel impulse re- 
sponse (CIR) h; = [h'^{0), /I'Jl), . . . , /I'JL - 1)]'^ of length 
L (in sampling periods). We denote by 9k the timing error of 
the A;th RSS expressed in sampling intervals Ts, while e^ is 
the frequency offset normalized to the subcarrier spacing. As 
discussed in (E], during IR the CFOs are only due to Doppler 
shifts and/or to estimation errors and, in consequence, they 
are assumed to lie within a small fraction of the subcarrier 
spacing. Timing offsets depend on the distance of the RSSs 
from the BS and their maximum value is thus limited to the 
round trip delay from the cell boundary. In order to elimi- 
nate interblock interference (IBI), we assume that during the 
ranging process the CP length comprises Nq > ^'^ax + L 
sampling periods, where 0max is the maximum expected tim- 
ing error. This assumption is not restrictive as initialization 
blocks are usually preceded by long CPs in many standard- 
ized OFDMA systems. 

2.2. System model 

We denote by X„(g) = [X„(ig), X™(z, + 1), . . . , X„(i, + 
V — 1)]^ the discrete Fourier transform (DFT) outputs cor- 
responding to the qth tile in the mth OFDMA block. Since 
DSSs have successfully completed their IR processes, they 
are perfectly aligned to the BS references and their signals do 
not contribute to X„j(g). In contrast, the presence of uncom- 
pensated CFOs destroys orthogonality among ranging signals 
and gives rise to interchannel interference (ICI). The latter 
results in a disturbance term plus an attenuation of the use- 
ful signal component. To simplify the analysis, in the ensu- 



ing discussion the disturbance term is treated as a zero-mean 
Gaussian random variable while the signal attenuation is con- 
sidered as part of the channel impulse response. Under the 
above assumptions, we may write 



K 



X„^i^q+v) = V [C I ^e^"''^"^^ Hk{ek,ek,iq+v)+Wrn{iq+v) 



fe=i 



(2) 

where ujk — 27Tek/N, Nt ^ N + Nq denotes the duration 
of the cyclically extended block and Cg^ is the code matrix 
selected by the A:th RSS. The quantity Hk{0k,£k,n) is the kth 
equivalent channel frequency response over the nth subcarrier 
and is given by 



Hk{ek.ek,n) = jN{ek)H'k{n)e-^^-^'>'/^ 



(3) 



where i/[.(n) = Efco^U^)^"^^''"^^^ is the true channel 
frequency response, while 



7Jv(e) 



sin(7r£) 
N sm{T:e/N) ' 



jT,e{N-l)/N 



(4) 



is the attenuation factor induced by the CFO. The last term 
in (|2]i accounts for background noise plus interference and is 
modeled as a circularly symmetric complex Gaussian random 
variable with zero-mean and variance CT^ — af^+aj^j, where 
a^^ and a-jci ^"^ '■^^ average noise and ICI powers, respec- 
tively. From ^ we see that 9k appears only as a phase shift 
across the DFT outputs. The reason is that the CP duration is 
longer than the maximum expected propagation delay. 

To proceed further, we assume that the tile width is much 
smaller than the channel coherence bandwidth. In this case, 
the channel response is nearly flat over each tile and we may 
reasonably replace the quantities {H'f.{iq + i')}^Jo^ with an 
average frequency response 



v-i 



^fc(9) = -^E^^(*^+«) 



V 



(5) 



Substituting ^ and (O into (|2|i and bearing in mind (|5]l, yields 



K 



Xr,S, +v)=Y, e^2-("«'=+-'^)5fe(<z) + w„,{iq + v) (6) 



k=l 



where Sk{q) — ^N{£k)Hf.iq)e J^'^^i^k/N ^jj^ ^g have de- 
fined the quantities 

4 , SkNT 



and 



rjk 



M -1 N 

4 ek 



(7) 



(8) 



V -I N 

which are referred to as the effective CFOs and timing errors, 
respectively. 

In the following sections we show how the DFT outputs 
{Xmiiq + v)} can be exploited to identify the active codes 
and to estimate the corresponding timing errors and CFOs. 



3. ESPRIT-BASED ESTIMATION 

3.1. Determination of the number of active codes 

The first problem to solve is to determine the number K of 
active codes over the considered ranging subchannel. For 
this purpose, we collect the [iq + v)t\\ DFT outputs across 
all ranging blocks into an M-dimensional vector Y(ig „) = 
[Xoiiq + v),Xi{iq + u), . . .,XM-i{iq + v)^ given by 



K 

E 

fc=i 



Y(«g,.) = VeJ'^'""=5fc(q)eM(a) + w(z,,„) (9) 



where vif(ig,„) = [wQ{iq+v),wi{iq+v), . . . ,WM-i{T'q+v)V 
is Gaussian distributed with zero mean and covariance matrix 
allM while eM(0 = [1, e^'^"?, e^'4-«, . . . , eJ"2-(Af-i)«]r. 

From the above equation, we observe that Y(ig,j) has 
the same structure as measurements of multiple uncorrected 
sources from an array of sensors. Hence, an estimate of K can 
be obtained by performing an eigendecomposition (EVD) of 
the correlation matrix Ry = E{Y(iq ,u)Y^(ig,„)}. In prac- 
tice, however, Ry is not available at the receiver and must 
be replaced by some suitable estimate. One popular strategy 
to get an estimate of Ry is based on the forward-backward 
(FB) principle. Following this approach, Ry is replaced by 
Ry = ^(Ry + JRyJ), where Ry is the sample correlation 
matrix 

^^ = nT7 E E Y(z,^„)Y^(z,,.) (10) 

while J is the exchange matrix with I's on the anti-diagonal 
and O's elsewhere. Arranging the eigenvalues Ai > A2 > 
■ • • > Am of Ry in non-increasing order, we can find an es- 
timate K of the number of active codes by applying the MDL 
information-theoretic criterion. This amounts to looking for 
the minimum of the following objective function ||6l 

j:{k) = ]rk{2V - k) \n{MQ) - MQ{V - k) \np{k) 

(11) 
where p{K) is the ratio between the geometric and arithmetic 

means of {A^+p A^f+a' • ■ ■ ' -^a^}- 

3.2. Frequency estimation 

For simplicity, we assume that the number of active codes has 
been perfectly estimated. An estimate of ^= [^1, ^2, ■ • ■ , £,k\ 
can be found by applying the ESPRIT algorithm to the model 
(|9]l. To elaborate on this, we arrange the eigenvectors of Ry 
associated to the K largest eigenvalues Ai > A2 > • • • > \k 
into an M X iiT matrix Z = [zi Z2 • • • z/^]. Next, we consider 
the matrices Zi and Z2 that are obtained by collecting the first 
M ~ 1 rows and the last M — 1 rows of Z, respectively. The 
entries of ^ are finally estimated in a decoupled fashion as 



ZTT 



1,2, 



,K 



(12) 



where {py{\), Py{2), . . . , py{K)} are the eigenvalues of 

Zy - (Zf Zi)-iZf Z2 (13) 

and arg{yOj,(fe)} denotes the phase angle of Py{k) taking val- 
ues in the interval [— tt, tt). 

After computing estimates of the effective CFOs through 
(fT2] i. the problem arises of matching each £_k to the corre- 
sponding code C^^ . This amounts to finding an estimate of 
ik starting from ^/j. For this purpose, we denote by |£max| the 
magnitude of the maximum expected CFO and observe from 
^ that (Af - l)^k belongs to the interval Ii,^ = [ik -P;ik + 
13], with 13 == le^axl Nt{M - l)/iV. It follows that the effec- 
tive CFOs can be univocally mapped to their corresponding 
codes as long as /3 < 1/2 since only in that case intervals 
{I(^ }fe=i ^e disjoint. The acquisition range of the frequency 
estimator is thus limited to |emax| < N/{2Nt{M — 1)) and 
an estimate of the pair {ik,£k) is computed as 



4 = round (M - l)e 



and 



£k 




h. 



M 



(14) 



(15) 



It is worth noting that the arg{-} function in ( fT2b has an in- 
herent ambiguity of multiples of 27r, which translates into a 
corresponding ambiguity of the quantity ik by multiples of 
M-1. Hence, recalling that Ik £ {0, 1, ... , K„^ax - 1} with 
^max < M, a refined estimate of £k can be found as 



oiF) _ 



[4: 



M-l 



(16) 



where [a;]j\/-i is the value of a; reduced to the interval [0, M — 
2]. In the sequel, we refer to ( fTSl ) as the ESPRIT-based fre- 
quency estimator (EFE). 



3.3. Timing estimation 

We call X„,(q) = [X„,{ig), X,n{iq + 1), • • • , X,n{iq + V- 
1)]'^ the F-dimensional vector of the DFT outputs correspond- 
ing to the ^th tile in the mth OFDMA block. Then, from ^ 
we have 



X™((j) 



K 

E' 
fc=i 



;^"2-™«'=5fe(g)ev(r7fc)+w„(g) 



(17) 



Wherew„j(g) = [Wm{iq),Wjrtiiq + l), ■ ■ ■ ,Wm.{iq + V -l)]'^ 

is Gaussian distributed with zero mean and covariance ma- 



j27ri) pjiTTV 



„j2n{V-l)n]T 



ti-ix cr^Iy while ey (ry) = [1, e 
Since X^ (q) is a superposition of complex sinusoidal signals 
with random amplitudes embedded in white Gaussian noise, 
an estimate of r]= [ryi, 7^2, ■ • • , Vk] can still be obtained by 
resorting to the ESPRIT algorithm. Following the previous 
steps, we first compute Rx — 5 (Rx 



JR^ J) with 



M-l( 



R 



X = 



MQ 



^ ^X™(g)X^(g). (18) 



m=0 q=0 



Then, we define aV x K matrix U = [ui U2 • • • u^^:] whose 
columns are the eigenvectors of Rj^ associated to the K largest 
eigenvalues. The effective timing errors are eventually esti- 
mated as 



flk = ^a-rg{Px{k)}, A; =1,2,. 



,K 



(19) 



where {px{l),Px{2), . . . , Px{K)} are the eigenvalues of 

Vx = (Uf Ui)-iUf U2 (20) 

while the matrices Ui and U2 are obtained by collecting the 
first V — 1 rows and the last V — 1 rows of U, respectively. 

The quantities {'fik}k=i ^^ eventually used to find esti- 
mates {£k, Ok) of the associated ranging code and timing er- 
ror. To accomplish this task, we let a — 0max(l^ ~ l)/(2^)- 
Then, recalling that < 6'fe < ^max, from ^ we see that 
{V — l)'i]k + a falls into the range 7^,, = [£h — a;£k + a]- If 
6'max < N/{V — 1), the quantity a is smaller than 1/2 and, 
in consequence, intervals {I(^ }t=i ^6 disjoint. In this case, 
there is only one pair {£k, 0^) that results into a given rjk and 
an estimate of {£k^ Ok) is found as 



£k =round((F-l)7}fc-^a) 



and 



N 



V-1 



Vk 



(21) 



(22) 



As done in Sect. 3.2, a refined estimate of £k is obtained in 
the form 



^Y^ 



Kfc V-l- 



(23) 



In the sequel, we refer to (|22] | as the ESPRIT-based timing 
estimator (ETE). 



3.4. Code detection 



From (\M and ( |23] l we see that two distinct estimates l^. ' and 



fif) 



are available at the receiver for each code index £k. These 



estimates are now used to decide which codes are actually ac- 
tive in the considered ranging subchannel. For this purpose, 

we define two sets I^f) = {ik^}k=i and /(-^^ = {4^^}f=i 
and observe that, in the absence of any detection error, it 
should be I^^^^ — /^^^ — {£k}^^i- Hence, for any code 
matrix Cm we suggest the following detection strategy 



ifme/(/)n/(^' 
if m i i^f^ n /(-^^ = 



4> Cm is declared detected 
Cm is declared undetected. 



(24) 



In the downlink response message, the BS will indicate only 
the detected codes while undetected RSSs must restart their 
ranging process. In the sequel, we refer to (|24] | as the ESPRIT- 
based code detector (ECD). 



4. NUMERICAL RESULTS 

The investigated system has a total of A^ = 1024 subcar- 
riers over an uplink bandwidth of 3 MHz. The sampling 
period is Ts = 0.33 ps, corresponding to a subcarrier dis- 
tance of l/{NTs) = 2960 Hz. We assume that R^ A sub- 
channels are available for IR. Each subchannel is divided into 
Q — IQ tiles uniformly spaced over the signal spectrum at 
a distance of N/Q = 64 subcarriers. The number of sub- 
carriers in any tile is V — A while M — 4. The discrete- 
time CIRs have L — \2 taps which are modeled as inde- 
pendent and circularly symmetric Gaussian random variables 
with zero means and an exponential power delay profile, i.e., 
E{\hk{i)\^} = crl •exp(-£/12)for^ = 0,l,..., 11, where 
af^ is chosen such that £'{||hfc|| } = 1. Channels of differ- 
ent users are statistically independent of each other and are 
kept fixed over an entire time-slot. We consider a maximum 
propagation delay of 6'max — 204 sampling periods. Rang- 
ing blocks are preceded by a CP of length Nq — 256. The 
normalized CEOs are uniformly distributed over the interval 
[—17, ii\ and vary at each run. Recalling that the estimation 
range of EFE is \ek\ < N/{2NTiM - l)),v^esetn < 0.1. 



10" W 













ECD 




















5;? 


\ i:::::::::1:::::::::t 


:::::::t:::::::::4:::;:.:...J... 
i I 


K=3 



/n-' H- i ^4 '4 ' i i f i i i \ i 



10 n 14 16 IS 20 22 



Fig. 1. Pf vs. SNR for iC = 3 when fl is 0.05 or 0.1. 

We begin by investigating the performance of ECD in 
terms of probability of making an incorrect detection, say Pf. 
Fig. 1 illustrates Pf as a function of SNR = 1/cr^. The 
number of active RSSs is 3 while the maximum CFO is either 
fl = 0.05 or 0.1. Comparisons are made with the ranging 
scheme proposed by Fu, Li and Minn (ELM) in Q. The re- 
sults of Fig. 1 indicate that ECD performs remarkably better 
than FLM. 

Fig. 2 illustrates the root mean-square error (RMSE) of 
the frequency estimates obtained with EFE vs. SNR for K = 
2 or 3 and il — 0.05 or 0.1. We see that the accuracy of EFE 
is satisfactory at SNR values of practical interest. Moreover, 
EFE has virtually the same performance as K or fi increase. 

The performance of the timing estimators is measured in 
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Fig. 2. RMSE vs. SNR for if = 2 or 3 when fl is 0.05 or 0.1. 
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Fig. 3. P{e) vs. SNR for if = 2 or 3 when il is 0.05 or 0.1. 

terms of probability of making a timing error, say P{e), as 
defined in |8 1. An error event is declared to occur whenever 
the estimate 9k gives rise to IBl during the data section of 
the frame. This is tantamount to saying that the timing error 
Sk — Sk + {—Nc,D + L)/2 is larger than zero or smaller than 
—Nq.d + L — 1, where Ng,d is the CP length during the data 
transmission phase. In the sequel, we set Ncd = 32. Fig. 3 
illustrates P{e) vs. SNR as obtained with ETE and ELM. The 
number of active codes is if = 2 or 3 while 51 — 0.05 or 0.1. 
We see that ETE provides much better results than ELM. 

5. CONCLUSIONS 



Compared to previous works, the proposed scheme exhibits 
increased robustness against residual frequency errors and can 
cope with situations where the CEOs are as large as 10% of 
the subcarrier spacing. 
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We have presented a new ranging method for OEDMA sys- 
tems in which uplink signals arriving at the BS are impaired 
by frequency errors in addition to timing misalignments. The 
synchronization parameters of all ranging users are estimated 
with affordable complexity through an ESPRJT-based approach. 



